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Abstract 

After discussing the intrinsic ambiguity in determining the hght quark 
mass ratio mu/frid, we reexamine the recent proposal that this ambiguity 
can be resolved by applying the QCD multipole expansion for the heavy 
quarkonium decays. It is observed that, due to instanton effects, some 
matrix elements which have been ignored in previous works can give a 
significant contribution to the decay amplitudes, which results in a large 
uncertainty in the value of rriu/md deduced from quarkonium phenomenol- 
ogy. This uncertainty can be resolved only by a QCD calculation of some 
second order coefficients in the chiral expansion of the decay amplitudes. 



It has been observed by a number of authors [1 — 5] that second order corrections 
in chiral perturbation theory can significantly affect the estimate of the hght quark 
masses. In particular, it was pointed out that the determination of irtu/md suffers 
from a large uncertainty due to the instanton- induced mass renormalization [^ Q: 

M -^ M{uj) =M + ujMi (1) 

where the real matrix M = dmg{mu,md,ms) denotes the light quark masses in 
the QCD lagrangian, M/ = ^(detM^)(M"'")~^ = ■^[rridrria^rriurria^rriurrij) is the 
instanton-induced second order mass with the pion decay constant / ^ 93 MeV, and 
L(j is a dimensionless parameter of order unity. Most of the previous analyses on M 
do not distinguish M from M , and thus the corresponding results can be interpreted 
as those on M for an arbitrary value of tu, which leads to a large uncertainty in the 
extracted value of mu/md. 

Recently it was argued that the above mentioned difficulty can be overcome by 
noting that the instanton-induced mass Mj is distinguished from the bare mass M 
through its 6'-dependence where 6 denotes the CP violating QCD vacuum angle [0. If 
one keeps the 6'-dependence explicitly, Mj always appears with the phase e*^ due to the 
winding number of instantons. Note that M and e^^Mj have the same transformation 
under the QCD chiral symmetry SU{3)l x SU{3)ii x U{1)a under which 

M -^ e'^LMR^ e^e + 3a, (2) 

where L E SU{3)l, R G SU{3)r, and a generates the anomalous U{1)a rotation. 
Then one may be able to measure M directly, not M involving an arbitrary unknown 
parameter u, by probing the ^-dependence of the QCD dynamics. Among quantities 
that probe the 6'-dependence, the matrix elements ((f)\GG\0) (0 = 7r° or rj) were 



considered in ref. [Bl. It was then argued that the ratio Ra = i .^^. / can be 

{ri\GG\0^ 

reliably determined by applying the QCD multipole expansion for the quarkonium 
decays ip' -^ J /ip + 0, which would allow us to precisely determine mu/md. In 



this paper, we wish to reexamine whether the QCD multipole expansion apphed 
for the quarkonium decays can provide a way to determine mu/md without doing 
any nonperturbative QCD calculation. Our result then confirms the conclusion of 
ref. 0, viz. in order to precisely determine mu/irid, one needs to calculate a QCD 
matrix element which receives a potentially large instanton contribution. Since such 
a calculation is not available at present, a rather wide range of rriu/ma, including 
zero, should be considered as being consistent with our present knowledge of QCD. 

To proceed, let us briefly review the points that were discussed in refs. |^, |^, 0]. 
In order to extract information on M, one usually considers measurable quantities 
whose M-dependence can be deduced from an effective lagrangian of hadron fields. 
This is not a severe limitation since it is hard to imagine a measurable quantity 
which may provide useful information on M through its M-dependence but can not 
be described by any hadronic effective lagrangian. To satisfy the Ward identities of 
the QCD chiral symmetry, the effective lagrangian is required to be invariant under 
the chiral transformations of the involved hadrons and also those of M and 6 given 
in eq. (2). Since M and e*^M/ have the same chiral transformation property, for 
any term in the effective lagrangian which is first order in M, e.g. Oi = ai (MQi), 
there exists also the second order term O2 = bic'^^ (MiQi) , where (Z) = tr(Z), a^ 
and bi are chiral coefficients which are calculable within QCD, and Qi is a generic 
local functional of hadron fields which is transformed under the chiral symmetry as 
M~^. This can be understood within QCD by noting that e^^Mj corresponds to the 
effective current mass induced by instantons [0, ^. For any QCD diagram which 
contains an insertion of M and thus gives a contribution to Oi, one can replace M 
by the instanton-induced effective mass e^^Mj to make the new diagram contribute 
to O2. This means that the M-dependent part of the effective lagrangian can always 
be written as 



where the elhpsis denotes other possible higher order terms. 

To be invariant under the anomalous U{1)a symmetry, the effective lagrangian is 
^-dependent in general. One may then expand the 6'-dependent effective lagrangian 
around ^ = 0: 

n=oo 

^effie) = E ^"^&/- (4) 

n=0 

Clearly C^jjr = Ceff{0 = 0) describes the normal CP conserving strong interactions 
while the terms of nonzero n describe the ^-dependence of the QCD dynamics, includ- 
ing the CP violation due to a nonzero 6. In view of the arguments leading to eq. (3), 
for any C],rK that includes the corrections of 0{M'^), one can define a transformation 
of the form: 

M ^ M = M + uMi, bi^bi- uai, (5) 

under which C^^rK is invariant^ The above transformation mixes the instanton- 
induced mass Mj with the bare mass M, and thus will be called as the instanton 
transformation (IT) in the following discussion. Obviously Cjlj {n ^ 0) in eq. (4) 
are not invariant under the IT. This is nothing but to mean that -Cg/j does not distin- 
guish the instanton-induced mass from the bare one, while Cjlj- (n ^ 0) distinguish 
since they arise from the ^-dependence. Based on this observation, it was stated in 
ref. PI that the normal CP conserving strong interactions which are described by 
£g/j- always concern the effective mass M with an unspecified value of cu, while the 
CP violating amplitudes or the nonderivative axion couplings probe directly the bare 
mass M. (Note that up to a small mixing with mesons, 6 corresponds to the con- 
stant mode of axion in axion models.) For instance, if det(M) = 0, then 6 can be 
rotated away (even in the case without axion) and axion becomes massless although 
det(M) ^ 0. 

Although jCj^f/s {n ^ 0) in eq. (4) might be useful also, one usually considers only 



^ Throughout this paper, whenever we say about the IT, it is assumed that the corrections of 
O^M"^) are included, while those of 0(Af'^) are ignored. 



>Cg A in extracting information on M from experimental data. Of course the reason is 
that it is quite nontrivial to find a hnk between the available experimental data and 
the terms of nonzero n. The IT was defined as acting on M and the second order 
chiral coefficients fej's. Since -Ce// is invariant, all measurable quantities described by 
£gjp i.e. expressed in terms of M and other parameters that appear in Ceff, are 
also invariant under the IT. This can be understood also by observing that the IT 
can be considered as a kind of renormalization group (RG) transformation associated 
with the instanton-induced mass renormalization [j^. In taking this specific instanton 
effects into account, one can use the following scheme. For 6 which is fixed to be 
zero, the contribution to C),rr from instantons with size p < pj^ is taken into account 
by redefining the chiral expansion parameter as M — ;> M + uMj, while that from 
larger instantons (p > pj^) appears in the second order chiral coefficients 6j. In this 
scheme, the transformation parameter u of the IT, being naturally of order unity, can 
be identified as uj{pi) = /^ ^ D{p), where -D(/i) denotes an appropriately normalized 
dimensionless instanton density. Then the IT corresponds to a renormalization group 
transformation^ changing p/, and all measurable quantities deduced from C].jj should 
be independent of our choice of fii, i.e. are invariant under the IT. 

As was argued previously [^, if we restrict ourselves to jC-efj, the invariance of 
Cl,rK under the IT prevents us from precisely determining rriu/md. All equations 
deduced from £gyj are covariant under the IT. As a result, any quantity which is 
sensitive to the IT, i.e. its variation under the IT for u = 0(1) can be comparable to 
the expected central value, can not be fixed by the invariant experimental data. (Of 
course the invariant combination of such quantities can be fixed.) One can not avoid 
an uncertainty whose size is characterized by the variation under the IT. For iriu/md, 



^ In considering the IT as a RG transformation, it must be noted that the scale /i/ is introduced 
only for the instanton-induced mass renormalization, but not for other kind of instanton effects. 
This isolation of the instanton-induced mass renormalization from other instanton effects can be 
easily achieved in the dilute instanton gas approximation. Then the conclusion that measurable 
quantities deduced from £j, A are invariant under the IT can be understood by restricting ourselves 
to the effects of relatively small instantons for which the dilute gas approximation is valid. 



its variation under the IT is roughly ums/Airf which can be as large as about 1/2, 
implying a rather large uncertainty. 

To be more specific, let us consider the case of using the pseudoscalar meson 
masses as experimental input. The relevant part of the effective lagrangian is P] 

+ Lr {x^U - xU^y + Ls {x^Ux^U + xU^xU^) , (6) 

where x = '^BqM for Bq = — (0|-uu|0) / p defined at chiral limit, and the SU{'i)- 
valued U denotes the pseudoscalar meson octet. For the above terms, the IT of eq. 
(5) can be written as P, ||: 

ruu ^ TUu-l — -^- (cyclic m u, d, s), 

AttJ 

Li-> Li -^^, (7) 

1287rfio 

where h^ = hi = 1 and hg = —2. Note that rus can be considered to be invariant 

since its variation is negligibly small compared to rris for u of order unity. In refs. 

g, g, it was found thatg 
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1}^±^ ^ {l + /XMfMl{Ml,-Ml, + Ml,-Ml,)-\ (8) 

rrid -mu 

where A^ = — 32L7i?ow^s//^ — c for c ~ 0.33 = Oinis/'i'Kf). These equations do not 

fix mu/md, but give only a curve on the plane of {mu/md, A) where A = irts/md for 

the first equation and \ = L-j for the second. Note that these curves are parametrized 

by tu, and the uncertainties in mu/md and A are given by their variations under the 

IT for uj = 0(1). We will encounter the same situation even when other kind of 

measurable quantities, e.g. the baryon masses and the decay amplitudes for r] — ^ 37r, 

^p' -^ J /ip + 7r°(77), are used in the context of an appropriate form of Ceff- 

•^ Throughout this paper, we will use nid/'ms ^ nis/^TTf, and thus ignore the corrections sup- 
pressed by either md/rris or md/4:Trf, while keeping only those of 0{ms/4:Trf). 



A definite value of mu/rrid would be obtained if one can choose a specific value 
of A, but it is possible only through a QCD calculation of A since A is sensitive to 
the IT and thus can not be fixed by experimental data alone. Such an attempt was 
made recently by Leutwyler |^ for \ = L-j. By invoking r^'-dominance, it was argued 
that L-j falls into a rather narrow range of negative values, thus ruling out rriu = 
in view of the second equation in (8). However it has been pointed out later that 
instantons significantly suppress the negative ?7'-contribution to Lj, while enhancing 
the positive contribution from the pseudoscalar octet resonances [|^. This would make 
r^'-dominance for Lj not valid, and thus results in a large uncertainty in the value of 
Lf, allowing rriu = 0. In fact, any quantity which is sensitive to the IT receives a 
potentially large contribution from instantons. This makes the QCD calculation of A 
and the precise determination of mu/md even more nontrivial. 

So far, we have argued that there exists an intrinsic ambiguity in determining 
^^u/'^d by usual manner. This ambiguity is due to the invariance of C),}f under 
the IT, and as was pointed out in ref. 0, might be resolved by including the non- 
invariant teims C/^L {n 7^ 0) in the analysis. In QCD, one can probe the ^-dependence 
through the matrix elements involving the insertion of the two gluon operator GG, 
e.g., U\GG\S) (0 = 7r° or r/) and (o|(GG')(GG')|o). The hadronic realizations of such 



matrix elements contain C}^)^ [n ^ 0) in general. In axion models, ( 0|GG|O ) describes 
the axion-meson mixing while the other matrix element is related to the axion mass. 
As a result, physical processes involving axion may provide information on these 
matrix elements, and thus on the bare quark mass M. However it is totally unclear 
whether the normal strong interaction data can also provide any information on these 
matrix elements. In ref. 0, it was observed that the matrix element ((l)\GG\Q 
appears in the QCD multipole expansion applied for the heavy quarkonium decays 
$' ^ $ + ($ = J/ip or T). The major purpose of this paper is to examine 
whether this allows us to determine niu/md without doing any nonperturbative QCD 
calculation. 



The quarkonium decays $' — i> $ + can be described by the effective lagrangian 
of eq. (3) with fij's containing the 5't/(3)-singlet heavy quarkonium fields $ and $' 
together with 0. The resuhing decay amphtudes H{(f)) (for $' at rest) can be written 
as: 

H[4>) = e,,fc$:$,p,. [x (M0) + 7 (M,0) + ...], (9) 

where $^ and $j denote the spin vectors of $' and $ respectively, pi is the momentum 
of = 0°A'^. Clearly the above amplitudes must be invariant under the IT of eq. (5), 
acting on M and also on 7 which is proportional to an appropriate combination of 6j's. 
As was discussed by Voloshin and Zakharov (VZ) 0, the QCD multipole expansion 
whose expansion parameter is the inverse of the heavy quark mass provides another 
expression for H{(j)): 

H{ct>) = e,,fc$:<l>,Xfc(0) W^o, (10) 

where Wq depends on the quarkonium wavefunctions, and 



X,(0) = (<P\g'ElDkB1\Q) ^ p,X{<p). (11) 

At the leading order in the multipole expansion, the quarkonium wavefunctions can 
be considered to be independent of M, and thus Wq is invariant under the IT. This 
means that X{(j)) which is written as 

X(0)=x(M0)+7(M,0) + ..., (12) 

is invariant under the IT of eq. (5), implying that 7 = 7/W0 is transformed as 

7—^7 — ujx. (13) 

In any case, one can define 

5,(0) ^ - U\g'BtD,E^,\0) = p,B{<P), (14) 



so that 

X(0)=A(0) + i?(0). (15) 

Note that A(0) can be written as A(0) = ^ (^I^^^Glo) where GG = Gl^G""^" . 
In the attempt to determine mujm^ using the quarkonium decays, the quantities of 
interests are Rx = ^k^'^M ^{^) and Ra = ^(7r°)/v4(?7). Then Rx can be fixed by 
the quarkonium decay data, while it is necessary to fix Rj^ to determine mu/frtd. In 
fact Ra = Rx at the leading order in M, however we need to include the corrections 
of 0{M'^) for a meaningful determination of mu/rrid. In ref. 0, following ref. |^, 
it was simply assumed that |-B(0)| ^ |y4(</))|, which would imply Rx — Ra even at 
OlM"^). Here we first argue that there is no reason for this assumption to be viable, 
and later show how the instanton-induced mass renormalization promotes Ra to be 
greater than Rx- 

The statement that B{(j)) may be significantly smaller than A{(f)) was first made 
by VZ [^ who observed that A{(f)) = 0{\/I^), and thus is enhanced by one power 
of Nc with respect to the naive A'"c-counting. This enhancement is due to the ?7o-pole 
where the 5'f/(3)-singlet r/o denotes rj' at chiral limit. Roughly we have 

{<P\g'GG\0) ^ {vo\9'GG\o) x Mj^„ x -^, (16) 

vo 

where the ?7o-pole (= 1/M^^) is 0{Nc). Note that the mass-squared mixing M^^^ 
between tjq and is suppressed by the hght quark mass M, but is 0(1) in the N^- 
counting, and (rio\g'^GG\0) = 0(l/-\/iVc) obeys the naive A^'c-counting rule. In fact, 
the same enhancement can occur also for -8(0) if {riQ\BfDkE^\0) is nonvanishing. 
Then -8(0) can be equally important as ^4(0) even in the large A'"c-limit. If -6(0) 
does not receive any contribution from the intermediate r^o and thus is 0(l/-\/iVc), we 
would have |-B(0)/y4(0)| ~ M3/A^, where A denotes a typical hadronic scale which is 
0(1) in the A"c-counting. Clearly then X(0) will be dominated by A{(f)) in the large 
Nc limit. However in the real case of Nc = 3, this is not necessarily true since M^/ is 
large enough to be comparable to A. Again -8(0) can be equally important as A{(j)). 



VZ noted also that the equation of motion DkE^ = igq^^q gives -8(0) an ad- 
ditional power of the QCD coupling constant g. However it is hard to imagine that 
this extra g means a real suppression of 5(0) compared to ^(0). First of all, g is 
essentially of order one, even for the renormalization point above ?7if,, although the 
loop suppression factor as/47r = g"^ /IQ-n"^ <^ 1. Note that the use of the equation of 
motion has nothing to do with the perturbative QCD loop expansion. Furthermore if 
we consider the valence quark contribution (in the sense of the parton model) to the 
vacuum to meson matrix element of an n-gluon operator, it would include a factor 
g'^. Before using the equation of motion, both v4(0) and -8(0) involve two-gluon op- 
erators. The equation of motion changes the two-gluon operator in -8(0) to B^q'^X°'q 
which includes only one gluon field. Then at least for the valence quark contribution, 
-8(0) is not higher order in g compared to ^(0) since the matrix element of B^q^X^q 
will be enhanced by g~^ compared to the matrix elements of two-gluon operators. 

We have argued that there is no a priori reason to expect that -8(0) is signif- 
icantly smaller than A(0). One must include -8(0) in the analysis, and then R^ 
can be significantly different from Rx at 0{M'^). Among the 0{M'^) corrections, the 
instanton-induced mass renormalization is of particular importance for the determina- 
tion of rriu/md- Thus from now on, we will concentrate on how the instanton-induced 
mass renormalization affects A{(f)) and B{(f)), so that promotes Ra to be greater than 
Rx- For this purpose, we study the IT of the second order chiral coefficients a and 
(3 that appear in the following chiral expansion^: 

v4(0) = a {M<j)) + a {Mjcj)) + ..., 

5(0)=6(M0)+/3(M,0) + ..., (17) 

where the ellipses denote other possible second order terms. Note that eqs. (12) 
and (15) imply that ^(0) and -8(0) can be expanded as above although these matrix 
elements can not be evaluated by using £gjj alone. 

'^ Here we ignore the electromagnetic effects which were shown to be neghgibly small p2| . 
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In order to derive the IT of a, let us express the physical meson fields 7r° and t] 



in terms of 03 and 
can be written as 



asr: vr 







rj = (j)s — e03. The mixing parameter e 



y/Snid-niu 



m. 



4vr/ 



1^ 



where k is a dimensionless parameter of 0(1) which represents the size of the 0{M'^)- 
corrections. Clearly e should be invariant under the IT since it describes the physical 
meson fields in terms of the original wavefunctions 03 and 08 which are untouched by 
the IT. This then gives the following IT of n: 



In fact, using the chiral lagrangian, e was evaluated as 

1287rSo 



(19) 



/ 



-(Lg + 3Lj) + chiral logs, 



(20) 



assuring the above IT of k (see the IT of Lj's in eq. (7)). 



Using eqs. (17) and (18), we can obtain 



^(^°) = o("^« -md)a 



47r/ 3 3a 



(21) 



while the anomalous Ward identity 



d^{qn^M-'q) 



J^tT{M-')G'''^^G% + 2qt^,q 



gives 



A(7r° 



427r^ rriu — rrid 



UMl 



0{ 



ATTf 



(22) 



3 uiu + md 

where /^ and M^ denote the decay constant and the mass of the physical vr'' respec- 
tively. Here and in what follows, ellipsis always denotes a dimensionless coefficient 
of 0(1) which is untouched by the IT. Then comparing (22) with (21) together with 

^In fact, due to the S'C/(3) breaking, the 7r°-(/)8 mixing can be different from the 77-^3 mixing at 
0{A'P). However this difference does not affect our analysis of the IT, and thus will be ignored. 
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the expressions of /^r and M^r given in ref . [^ , we find 



a 
a 
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(Lg + 3L7 + SLe) + 



(23) 



which yield the following IT 



a 



a + 2a;a, 



/3^/3-o^(6 + 3a). 



(24) 



Here the IT of (3 is derived from that of a and 7, using x = a + h and 7 = 0; + /?. 

We are now ready to see how the instanton-induced mass renormalization affects 
Ra and Rx-, and what is still necessary to determine rriu/md using the QCD multipole 
expansion applied for the quarknonium decays. Using eqs. (12), (17), and (18), we 
can obtain 



Rx 
Ra 



X( 



TX 



X{ri) 4 irts 

A{tx'^) 3\/3 md — rriu 



l + ^.ix 



1 + 



47r/ 
4nf 



U 



(25) 



A{ri) 4 

where ^x = '«/3 — 27/30; + ..., ^a = /^/3 — 2a/3a + ... are dimensionless coefficients 
of 0(1). Then the IT of a, 7, and n derived so far gives the following IT: 



^X ^ ^X+ UJ, U-^ ^A- ^■ 



(26) 



We have argued that the IT can be interpreted as a kind of RG transformation 
changing the scale /x/ that appears associated with the instanton-induced mass renor- 
malization. With this interpretation, the IT of S,a,x for a maximal value of the trans- 
formation parameter, i.e. u = cUmax, represents the negative of the full instanton 
contribution to C,a,x- Then C,a receives a positive contribution uj^^x from instantons, 
while ^x receives a negative one with the same size. Although its precise value is 
quite sensitive to the unknown nonperturbative QCD dynamics, it has been observed 
that Climax can be as large as 2 ~ 3 in the semiclassical instanton gas approximation 
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|T], ^. Then although not very rehable, for such a large cUmax, it is more likely that 
Ra is significantly greater than Rx- 

Applying the QCD multipole expansion for the measured quarkonium decay 
widths V{i)' ^ Jip + ^) (0 = 7r° or r/), one finds Rx ^ 4.3 x IQ-^ [g ^ ^. Also 
in ref. 0, the size of Ra = (^d + m„)_R^/(mrf — m„) was estimated within the sec- 
ond order chiral perturbation theory of the light pseudoscalar mesons, which yields 
Ra — 8 X 10~^. Note that both Rx and Ra are invariant under the IT, and thus 
can be fixed by experimental data. If we assume as in ref. |0 that |-B(0)| <^ 1^(0)1 
which implies Ra — Rx or equivalently C,a — C,x, the measured values of Rx and 
Ra would give niu/md ~ 0.3. However we already argued that there is no reason 
for 1-8(0)1 ^ |yl(0)|. Furthermore once we include -8(0) in the analysis as it must 
be, instantons give a potentially large positive contribution cUmax to ^a^ while ^x 
receives a negative contribution with the same size. This is essentially due to the 
instanton-induced 0{M'^) piece in B^-k^), i.e. the term with the coefficient (3 arising 
from the instanton-induced mass renormalization. Then what we obtain from the 
entire analysis can be summarized by the equation: 

' ^ 0.54 1 + ^ ^A - ix)). 27 

where 8^ = [^a ^ ^x) is a totally unknown coefficient of 0(1). If cUmax = 2 ~ 3 whose 
possibility was assured within the semiclassical instanton gas approximation |l|, ^ , it 
is conceivable to assume that 6C, is dominated by the positive instanton contribution 
2cUmax5 implying ^a > ^x- Then the measured values of Rx and Ra would give 
TTT'u/'tTT'd < 0.3. In any case, in order for mu/md to be precisely determined, we still 
need to compute (within QCD) the chiral coefficient 6^ which is very sensitive to the 
IT 

To conclude, we have examined whether the QCD multipole expansion applied 
for the quarkonium decays $' ^ $ can be useful for the precise determination 
of mu/md. The matrix element -8(0) which has been ignored in previous works 
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can significantly affect the estimate of mu/md, particularly through the 0(M^)-piece 
induced by instantons. As a result, a rather wide range of irtu/rnd (including zero) 
can be consistent with the observed quarkonium decays. As was concluded in ref. |^, 
to determine m^/m^ precisely, we need a QCD calculation of the chiral coefficient 
^^ = {^A ~ C,x) which is sensitive to the IT. If instanton contribution to 6C, dominates 
over other contribution, which is conceivable for cUmax = 2 ~ 3, and thus 6^ > 0, 
the observed quarkonium decays imply iriu/md < 0.3, allowing the massless up quark 
scenario [1, 2, 3] for the absence of CP violation in strong interactions. 
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